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Abstract 

This paper studies the chemotaxis-haptotaxis system 

{ ut = A u — xV • (uVo) — £V • (uVui) + pu(l — u — w), (a:, t) G fi x (0, T), 

v t = Av — v + u, (x,t) £ fl x (0,T), (*) 

wt = —vw, (x,t) € x (0 ,T) 

under Neumann boundary conditions. Here SI C I 3 is a bounded domain with smooth boundary 
and the parameters £, x, A 4 > 0- We prove that for nonnegative and suitably smooth initial data 
(uo,vo,wo), if x/M is sufficiently small, (★) possesses a global classical solution which is bounded in 
D x (0, oo). We underline that the result fully parallels the corresponding parabolic-elliptic-ODE 
system. 


1 Introduction 

The motion of cells moving towards the higher concentration of a chemical signal is called chemotaxis. 
A classical mathematical model for this type of processes was initiated by Keller and Segel in the 1970’s 
[5], and in a prototypical form this model writes 

f u t = Au -xV • (uVv), (i,i)eSlx(0,T), 

| Vt = Av — v + u, (x, t) £ fl x (0, T). 

Here u is the density of cells, v denotes the concentration of the chemical substance, and x > 0 measures 
the sensitivity of chemotactic response. In recent 40 years, a large quantity of literature has been devoted 
to study the global existence as well as singularity formation in either finite or infinite time for this system 
mmm- We also refer to ei Eiia for a broad overview. 

Apart from this classical Keller-Segel system, a large number of variants has been proposed to describe 
taxis phenomena in mathematical biology. Among them, a model for tumor invasion mechanism was 
introduced by Chaplain and Lolas [3j- In this model, tumor cells are assumed to produce a diffusive 
chemical substance, the so-called matrix-degrading enzyme (MDE), which decays non-diffusive static 
healthy tissue (ECM). It is observed that both the enzyme and the healthy tissue can attract the cancer 
cells in the sense that the cancer cells bias their movement along the gradients of the concentrations of 
both ECM and MDE, where the former of these processes, namely taxis toward a non-diffusible quantity, 
is usually referred as haptotaxis. 
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Additionally, the cancer cells compete for space with ECM, and at the considered time scales moreover 
logistic-type cell kinetics need to be taken into account. If futhermore the ability of ECM to spontaneously 
renew is included, the Chaplain-Lolas model becomes 

! ut = A u — \V • (uVw) — £V ■ (uV w) + fxu{ 1 — u — w), (x, t) £ Q x (0, T), 

rvt = Av — v + u, (i,t)e!lx(0,r), (1.2) 

wt = —vw + r/w(l — u — w), (x, t) £ fl x (0, T), 


where u, v and w denote the density of cells, the concentration of MDE and the density of ECM, 
respectively, where the parameters £, x, fx, ry are positive constants and r > 0, and where SI C 1^, 
N > 1, denotes the physical domain under consideration. 

Assuming w = o, Ol is reduced to the classical Keller-Segel system with logistic source, which has 
extensively been studied during the past 20 years. Compared with the pure chemotaxis system mentioned 
above, one may expect the logistic source and, especially, death terms to enhance the possibility of 
bounded solutions. In fact, Tello and Winkler [T5| proved that if r = 0 and 


M > 


(N - 2 ) + 

N 


X, 


(1.3) 


then for any regular initial data, the system admits a unique global classical solution which is bounded. 
In the case r = 1, it is known that bounded solutions exist in lower dimensions (N = 1, 2) for any fj, > 0 
@], and that the same result holds for ji > /jq with some Hq(x) > 0 in higher dimensions [20]. More 
precisely, a careful inspection of the proofs therein shows that in fact large values of the ratio are 
sufficient to exclude blow up in either finite time or infinite time. 

Concerning El with possibly nontrivial w, the strong coupling between remodeling and chemotaxis 
substantially complicates the situation, and accordingly the knowledge on this topic is quite incomplete 
so far. To the best of our knowledge, global existence of weak solution is obtained in m for N < 3, 
where El is included as a subsystem. And global solvability of classical solutions in this full system is 
known only when r = 0 and N = 2 m- Disregarding the chemotaxis effect, the haptotaxis-only version 
with x = 0, r = 1 was studied in [T2] . 

In the real situations, the ECM degrades much faster than it renews, thus the remodeling effect can be 
neglected, that is, we may assume rj = 0. Under this hypothesis, the corresponding parabolic-elliptic 
simplification r = 0 has been studied by Tao and Winkler in he where it has been proved that solutions 
stay bounded under the same condition as in the case w = 0, that is, when (11.311 holds. This shows that 
in this situation the haptotaxis term does not affect the boundedness of solution, and that accordingly 
the chemotaxis process essentially dominates the whole system. A natural question is whether a similar 
conclusion holds in the fully parabolic system obtained on letting r = 1. In he Tao gives a partially 
positive answer in this direction by proving that when N = 2, solutions remain bounded for any /x > 0, 
which thus parallels known results both for r = 0, and also for r = 1 when w = 0. As far as we can 
tell, however, despite a result on global existence established in [14] , the question of boundedness of 
solutions is completely open in higher dimensions. It is the purpose of this work to furthermore establish 
a corresponding parallel result for the three-dimensional parabolic-parabolic-ODE chemotaxis-haptotaxis 
model in this direction. 
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Accordingly, we deal with the system 


ut = A u — yV • (uVv) — ■ («Vw) + p,u( 1 — u — w ), 


vt = Ac — c + u, 

< Wt = —vw, 

= % = 0, 

u(x, 0) = uo(aj)) c(x, 0) = vq(x), w(x, 0) = wq(x), 


(x,t) t!!x (0,T), 
(x,f) g!1x (0,T), 
(x, l)ellx (0, T), 
(x,t) easix (0,T), 

X £ fl, 


(1.4) 


where fl C K 3 is bounded with smooth boundary and y. £, /z > 0. We assume that initial data are regular 
enough and satisfy a standard compatibility condition in the sense that 


u 0 £C 0 (fl), v 0 £ w 1 ’ 00 (ft), xco e c 2 ’“ (ii)(a e (0,1)), 
^ = 0. 


(1.5) 


Then our main result says the following. 


Theorem 1. There exists 9 q > 0 such that whenever y > 0, u > 0 and £ > 0 are such that ^ < 9 q, 

_ A 4 

for any initial data (uq, vq, u>o) fulfilling (11.51) . there exists a unique classical solution (u,v,w), which is 
global in time and bounded infix (0, oo). 

Remark 1.1. We only carry out the proofs for three-dimensional case in this paper. Actually, we can 

2 N 

start the iteration from the higher regularity of v £ L N ~ 2 (fl), the same result still holds for N < 15. 

We see that although our hypothesis on the parameters is not as explicit as m obtained for the 
parabolic-elliptic counterpart, it still shows that again boundedness of solutions is enforced by a condition 
merely referring to the interplay between chemotaxis and quadratic degradation in logistic source. 

Apart from this, we find it worth mentioning that our approach even shows a new result for the pure 
fully parabolic chemotaxis system with logistic source in the sense that when w = 0, N > 3, the system 
admits a classical bounded solution if ^ is sufficiently large. Compared with a similar conclusion under 
the alternative assumption that Ar be large 1201 . our result seems more consistent with m for the 
parabolic-elliptic system where the linear ratio ^ is found to determine the boundedness of solution. 


2 Preliminaries 

We first state a result on local existence of classical solutions. Without essential difficulties, the proof 
can be derived based on that in HU Lemma 2.1]. 

Lemma 2.1. Let N = 3, y > 0, £ > 0 and /z > 0. For (uq, Vo, wo) satisfying 11.5 1) . there is T max £ (0, oo] 
such that admits a unique classical solution 

u £ C°(ti x [0, T max )) n C 2,1 (D x (0, T max )), 

V £ C°(n x [0, T max )) n C 2,1 (fl x (0, T max )) n L“ c ([0, T max ); w 1 ’ q (fl))(q > 3), 

W £ C 2 ’ 1 ^ x [0, T max )), 

such that 

u> 0, u > 0 and 0 < w; < ||rco||i=°(n) for all te[0,T max ). (2.1) 
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Moreover, if T max < oo, then 


limsup(||u(-, t)||z,«>(n)) = oo. 

t/ A 'I max 


( 2 . 2 ) 


According to the above existence theory, we know that if we fix any so £ (0,T max ), then there exists 
M > 0 such that 


s o)IU“(n) + ll u (’> s o)||L°°(n) + \\ w ('j ■ s o)llw 2 ’°°(n) < M. (2.3) 

Observing that w can be represented by v and w(x, So), we can compute Aw in a convenient way. Upon 
a slight adaptation of m Lemma 2.2], we can prove a one-sided pointwise estimate for Aw as follows. 

Lemma 2.2. Let (uo,vo,wq) satisfy (11.51) . ( u,v,w ) solve (11.41) . For all So > 0, we have 

A w(x,t) > Aw(x,so) ■ e~^» v ^ a! ’ s ^ ds — 2e~^ s o v ^ x,s ^ ds Vw(x,so) ■ / Vv(x,s)ds 

Js 0 

~—w(x, so) — w(x, so)v(x, t)e~ ^ s o v ( x ’ s ^ ds (2.4) 

e 

for all x £ U and all t £ (so,T max ). 

Proof. Representing w(x,t) according to 

w = e s o w(x,so) (4-5) 

for all x £ Q and t £ (sq, T max ), we directly compute 


Aw(x,t) = Aw(x,so)e ^o v( - x,s ^ ds — 2e v< - x ’ s - lds Vw(x, so) • I' S7v(x,s)ds 

J s 0 

+w(x, so)e~ f =0 ' u ( 2 -’ s)d;s | f Vw(x, s)ds| 2 — w(x, so)e _ v ^ x ’ s ' )ds f Av(x, s)ds. 

j Sn J Sn 


Since ze z < by dropping some nonnegative terms, we obtain that 


Aw(x,t) > Aw(x,so)e J°o l ’ (x,s ^ ds — 2e v ^ x,s ' >ds Vw(x, so) • f 'Vv(x,s)ds 

•J So 

t r ^ 

— w(x, so)e~ ^ s ° v ( x ’ s ) ds / (v 8 (x, s) + v(x, s) — u(x, s)) 

j So 

> Aw(x,so)e~^ s o v ^ x ’ s ‘ >ds —2e~^‘o v ^ x ’ s - >ds Vw(x,so)- f Vv(x,s)ds 

J So 

— w(x, so)e _ ^ 3 ° V< ' x ’ s ' >ds (v(x,t) — v(x, so)) — w(x, so)e~ v ( x ’ s ) ds f ^(^ s )^ s 

Jo 

> Aw(x, So)e _ v(x,s)ds _ 2 e ~ L 0 v ( x ’ 6 ' >d6 Vw(x, so) • l Vv(x,s)ds 

Jo 

)(x, sq)v(x, t)e~ ^ s o v ( x < s ) ds — - w (x, sq) 


— wl 


for all t £ (so,T max ). Thus the proof is complete. 


□ 


With the aid of Lemma 12.21 we can furthermore prepare a preliminary estimate of an integral related 
to the haptotactic interaction. This estimate will be used in different ways later on. 
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Lemma 2.3. Let \ > 0, £ > 0, and assume that ESI) holds. Then for any p > 1, Sq G (0,T max ), the 
solution of satisfies 

(p-l)£ [ u p_1 Vti • Vm < (3M£ + -M£) [ u p + M£ [ u p v + 2M{p-l)f [ u^lViil (2.6) 
J i ^ J q, J n J q 

for all t G (s 0 ,T max ). 

Proof. Integration by parts and an application of Lemma 12.21 yield that 

(p — 1)£ f u p_1 Vzt • Vw 
Jn 

-£ f u p Aw 

J n 

-£ f u p (Aw(x, so)e~ f s o v ( x,s ) ds — 2e~ f‘ , o v ( x,s - >ds Vw(x, so) ■ f Vv(x,s)ds 

J f2 J so 


< - 


P ' Jn 
P~ !, 


P do 
1 


-wj(x, so) — iy(x, so)u(x, t)e ^o v ( x > s ) ds ) 

e 

< + [ u p + Mf f u p v-2?—^-C [ u p \7w(x,s 0 )-Ve- f ‘o v ( x ’ s)ds 

e do do P do 

= (Af£+-AfO [ u p + Mf [ u p v + 2^-^-f [ u p Aw(x,s 0 )e- f ‘o v(x ’ s)ds 
e do do P do 

1 0?3 —£ / Vm p -V w(x,so)e“ / “o ,,(a: ’ s)ds 
P do 


□ 


< (3M£ + -M£) [ u p + Mf f u p v + 2M(p-l)£ [ u^lVul 

^ Jn Jn Jn 

for all i G (so,dmax)) where in accordance with (12.31) . M is an upper bound for ||w(-, so)||vk 2 '°°( 0 )- 
Lemma 2.4. Let \ > 0, f > 0 and p > 0, and assume EH- Then there exists C(p, |fl|) > 0 such that 

[ u(;t)<C, [ V (;t)<C (2.7) 

do do 

for all t G (0, T max ). 

Proof. The first inequality in (11.41) can be proved by simply integrating the first equation in (11.41) on Q 
and using that (J n u) 2 < |f2|(J n u 2 ) due to the Cauchy-Schwarz inequality. The estimate of J n v can be 
obtained in a similar way and with the aid of the first inequality. □ 


As an essential ingredient of the proof of our main result, we will use a Maximal Sobolev regularity 
property associated with the second equation in (11.41) . The following lemma is not the original version of 
a corresponding statement in [5J Theorem 3.1], but by means of a simple transformation it is adapted to 
the current situation by including an exponential weight function. 

Lemma 2.5. Let r G (l,oo), T G (0,oo). Consider the following evolution equation 

v t = Av-v + f, (x,t) G Cl x (0,T), 

fp = 0, (x,t) € an x (0,T), (2.8) 

v(x,0) = vo(x), x G Q. 

For each vq G W 2 ' r (Ll) such that = 0 on dLl and any f G L r ((0, T); L r (Ll)), there exists a unique 
solution 

v G fT 1 ’ r ((0, T); L r (Ll)) n L r ((0, T); TT 2 ’ r (f2)), 
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and there exists C r > 0 such that 
r T 


[ f e^ rs \Av{x,s)\ r dxds < C r [ f e^ rs \f(x, s)\ r dxds + C r \\vo\\ r L r (Q) + C r \\ Ai;o||£ r(n) . (2.9) 

Jo Jn Jo Jn 

Moreover, if for some Sq £ (0 ,T), v(-,Sq) satisfies v(-,sq) £ W 2,r (fl) with = 0 on dLl, then with 

the same constant C > 0 as above we have 

f f e% rs \Av(x,s)\ r dxds < C r f j e? rs \f(x, s)\ r dxds + C r \\v(-, s 0 )\\ r L r {n) + C r \\Av(-, s o )\\ r Lr{n {2.10) 
J So J ^ J So J 

Proof. Letting w(x,s ) = e^ s v(x,s) — x( s ) v o, where x(s) is a cut-off function such that % £ Co°([0,1)), 
and x(0) = 1, we see that w solves 

w s (x,s) = (A - ±)w(x,s) + e? s f(x,s) +x{s)Av 0 - §x(s)t>o - x'( s ) w o, (x,s) £ flx (0 ,T), 


= 0, 


dw 
dv 

w(x , 0) = 0, 


(x,s) £ dn x [0,T), 

X £ ft. 


Let d be A - ) associated with Neumann boundary condition, we know that A is a generator of an 
analytic semigroup of negative exponential type. An application of standard results on maximal sobolcv 
regularity provides c r > 0 such that 

[ im-,s)[|£rds+ [ \\ws(-,t)\\ r L rds + f ||Au>(-, s)||£ r(n) ds 
Jo Jo Jo 

J o \\e^f(-,t)\\ r L r + ||x(s)Au 0 - ^x(s>o - x'(s)v 0 \\ r L rds 


< c 


< c r 


\e 2t f{-,t)\\ r L r + CrdWvoWw^.^n). 


( 2 . 11 ) 


where c\ := 2 + sup y'(.s). The triangle inequality with (12.111) implies that 

« 6 [ 0 , 1 ] 

J l|Au)(-,s)||^ds < j ||(A- ^)w{-,t)\\ r Lr ds + J ^\\w(-,t)\\ r Lr ds 

< c r f \\ei t f{-,s)\\irds + c 1 c r \MW’.r, 

Jo 

which leads to (12.91) upon direct computation and letting C r := max{c r , CiC r }. Thereafter, (12.101) directly 
follows upon replacing v(x, t) by v(x, t+s o), where the constant might be different from the aforementioned 
one. □ 


3 Boundedness 

In this section, we derive the claimed boundedness result via combining the above result on maximum 
Sobolev regularity with a Moser-type iteration. The former ingredient is first used to estimate u in some 
appropriate Lebesgue space, from which a certain suitable estimate of Vv will follow. This approach will 
be carried out twice to ensure that Vv is bounded in L°°(fl). Thereupon we can establish a series of 
inequalities based on which a Moser iteration is performed to finally achieve boundedness of u in L°°( fl). 

Lemma 3.1. Let X > 0, £ > 0, /i > 0, 9 = ^ and let p\ £ (|,2), pi £ (3, oo). Then there exists 9 q > 0 
such that whenever 9 < 9q, for any (u o, Vq,wo) fulfilling 11.51) and some Sq £ (0, T max ), we can find C > 0 
such that 

[ u P2 (-,t)<C for all t £ (s 0 ,T max ). (3.1) 

Jn 
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We are going to prove Lemma l3.1l bv several steps. Let us first provide an important ingredient for the 
estimate of ||zt(-, i) ||z,i»(si) - The next Lemma offers a general iteration step from v to u. We will use it to 
first estimate ||m(-, i)IU p i(n) with pi £ (|,2), and then ||u(-,i)llL p 2 (£ 2 ) with p 2 £ (3,oo). These results are 
under different conditions on —. For convenience, we assume that both of the conditions hold in Lemma 

A 4 

13.11 such that Lemma EZH is applicable for both pi and p 2 - 


Lemma 3.2. Let y > 0, £ > 0, p > 0, p > l. Then there exist constants 9 p > 0 and C > 0 such that if 
9 = — < 9 V and 

fi y 



for all t £ (s 0 ,Tmax)- 


Then we have 


t) < C for all t G (S 0 , Tmax)- 


(3.2) 


(3.3) 


Proof. First we see that for any a, b > 0, Young’s inequality provides k p > 0 such that 


ab < —a p + k p b p+1 . (3-4) 

8 

Let C r +i denote the constant from Lemma [2751 for r G (1, oo). Now we can find 9 P > 0 small enough such 
that 


Cp+ik p Q p+1 < — for all 8 < 6 P . 


Testing the first equation in (11.411 with u p 1 (p > 1) and integrating by part imply 

~ f vP + {p-T) [ u p ~ 2 \Vu\ 2 
P at J n J n 

= (jp — l)x [ u p ~ 1 \7u ■ Vv + (p — l)f f m p_ 1 Vm • Vw 
Jn Jn 


+ p u p — p / u p+1 — p u p w 
J Cl J Cl J Cl 


< 


p-1 


< - 


v 

p-1 


X f Vu p • Vv + (jp — 1)£ f u p 1 Vu-Vw + p f u p — p f 

J Cl J Cl J Cl J Cl 

X f u p Av + (p—l)f J u p ~ 1 Vu-\/w + p j u p — p f 

J Cl J Cl J Cl J Cl 


# p +1 


P J Cl J Cl J Cl J Cl 

for all t G (so, Tmax). We see that (12.61) and (13.171) entail the existence of C 3 > 0 satisfying 

(p — 1)£ f m p_1 Vu • Vtu 

Jn 

< c 3 f f u p + c 3 f f u p v + c 3 pf j u p_1 |Vw| 

J Cl J Cl J Cl 

< C3{ jf „ + g 1 + 1 + ^ J a + |^/ n » 

u p - 2 |Vu | 2 + k p p- p ^ +1 \\v(-,t)\\ p +l Hn) 


< (c 3 £ + 


)[* + £[ + 


n 

P~ 1 


2(P-1) 


for all t G (sojTmax). From (13.41) . we estimate for all t G (so,T ma x), 
p-l_ 


X [ u p Av < x [ w p |An| < ^ f u p+ 1 +k p x p+1 p p f |Au| p+1 . 
Jn Jn ° Jn Jn 


(3.5) 


(3.6) 


P Jn 


(3.7) 


(3.8) 
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Inserting (13.711 and (13.81) into (13.61) and some rearrangement yield 


1 d 
V dt 


P~ 1 


'n 




— “tM / u 


p+i 


4 

p+ 1 
~^P~ 


n 


c 3 C + +^)J n uP + KP~ P X v+l J n |A^| p+1 + k p p- p e +1 \W;t)\\ P L + P l 


(Q) 


P + 1 


C |C 2 P 2 
+ c 3 £ + 3S 


2(P-1) 


+ p I I u p --a \ u p+1 


3 

4' Jn 


+fc pX ^+v- p [ lA^r 1 + k P v- p e +1 \H-M p L + P l m 

j q 


for all t £ (sojTmax). We again apply Young’s inequality to obtain that 

'p + 1 , ^ 4( 2 p 2 


(Hr + + W=T) + ' ") L "" - 4 In “” + ‘ + 


2 p 1 ' 2 (p - 1 ) 

where C 5 > 0 depends on p,£,p. Upon (13.91) and (13.101) . we infer that 

'p + 1 


(3.9) 


(3.10) 


T / «”<- 
dt in 


77 P [ u p+1 + pk P x p+1 p p I |Au | p+1 
^ Jn Jn 


+ c 4 (/z, £,p)p + pk p p p £, p+1 ||u(-,t)||ilii (n) 

for all t £ (sojTmax). Applying Gronwall inequality in different form to the above inequality shows that 
f u p (-,t) < e -(4 i )0--o) f «*(., ao ) - /’ t e -(4 i )0-U f u p+1 {;s)ds 

J Q J Q " J so JQ 


+ kpx P+1 p p f e ^ 2 i 4 f |Au(-, s)| p+1 ds + f e ( P2 s )(c 4 (/z,£,p) 

« So ** ^ J So 

+ V P e + ik(-U)ll^ 1 ( n)) rfs 

= e -(4 i )(^o) r uP{ ^ so) _p-(z±±) t r [ e( 4 i ) v+1( , s) 

./ fi " J so J Q 

+ pkpX p+l p- p e-^y [ [ e( £ ^ i ) s |Az;(-,s)|P +1 
J So J ^2 

+ (c 4 (/z,^p)+fc p /z- p ^ +1 ||z;(-,t)||^i(n)) 


(3.11) 


for all t £ (so,T m ax)- In order to estimate the third term therein, let us note that an application of 
Lemma 12.51 results in 


pk P x P+1 p p e ( P 2 Y f j e ( P 2 ) s |Au | p+1 
J So " fi 

< pkpC p+1 x p+1 p~ p e~^y J^ E i i yu p+1 + C p+1 pkpX p+1 p~ p e~^y\\v(-,s 0 )\\ p ^l P+1{n) 

= pk p C p+ lX p+1 p- p e~^y f [ e^>u p+1 + Cp +1 pk P x p+1 p- p e-( s ^) t M p+1 (3.12) 

J Sq J fi 


for all t £ (sojTmax) and M as in (12.31) . Since 9 = combining 


I, we finally arrive at 


f u p (;t) < [ U p (;s 0 ) - pp (l ~ k p C p+1 e p+1 ) e-^> f [ e(^) 

Jn Jn / Js 0 Jn 

+ C p+1 pkp X p+1 p- p e-( s ^) t M p+1 


s u p+1 














+ (c4(M,p) + k^- p e +1 H-,t)\\ P L + P l ( n)) 

for all t G (so,T max ). We see from (13.21) and the condition on 9 p in (13.51) that 

/ u p (-,t) < C(n,x,£,P, lk(' ; So)||^2, P+ i (n) ) 

J n 

for all t G (so,T max ) upon obvious choice of C > 0. Thus the assertion is derived. 

Now we can apply Lemma [7T?1 to improve the regularity of u from L 1 ( ft) to L Pl (fl). 
Lemma 3.3. Let pi G (|,2) and 9 < 6 pi , there is a constant C > 0 such that 

||w(-,t)||iPi(n) < C for all te (s 0 5 ^max) • 

Proof. First we use the representation formula for u, 


(3.13) 

(3.14) 

□ 


j(-,t)=e tA v o+ f S ^ A 1 ' > u(-,s)ds 
Jo 


(3.15) 


(3.16) 


for all t G (0,T max ). Let p\ G (|,2), the L p — L q estimates for the Neumann heat semigroup [19} 


Lemmal.3] and (12.71) allow us to pick constants Ci > 0 and C 2 > 0 such that 

f || e ( *— s) ( A—1 )'ii(-, s 


IK 


TJIIlpi+^o) A 


< ||e t(A_1) 


LP l + !(fi) 


ds 


^01| L p i + ! (Q) + J 

< e *ll u o||i,pi+i(n) + ci f (1 + (t — s) p i+^)e ^ s - ) || m(-, s)Hii(n) 
Jo 


ds 


(3.17) 

□ 


< c 2 

for all t G (0,T max ). An application of Lemma I53il and (13.171) lead to (13.151) . 

With the higher regularity for v obtained in the above lemma, a similar reasoning for Lemma 13.21 will 
provide us higher regularity for u. Now we are ready to prove Lemma 13.II 

Proof of Lemm,a, \3.1\ Let p\ be choosen as in Lemma [3.21 and p 2 G (3,oo), define 0q = min{0 Pl , 0 P2 }. 
According to Lemma 13.21 and standard estimates for Neumann heat semigroup generated by A, we see 
that 

IK-,*)||z/~(n) < ||e t(A_ 1 ) u(-,s 0 )|| L =c ( o ) + / || e ( ‘-*)( A - 1) u (. ; s)\\ L ^ {Q) ds 

J 0 


<e i ||w(-,so)||i=o(o)+ci / (1 + (t-s) 2pi )e { * s) ||u(-> s)\\ LP1 ^)ds 


(3.18) 


for all t G (s 0 ,T max ) with some C\ > 0. Lemma [3.31 guarantees the boundedness of the above estimate. 
Also Holder’s inequality yields that 


H^(‘> t) IIL p 2 + 1 (f2) < Cl, t G (sq, 7 max ). 


(3.19) 


Due to the definition of do and (|3.19|) , now we can apply Lemma 13.21 to find the existence of a constant 
C(p,x,£,P 2 , IK>so)||w 2 ’ p 2 +i(n)) > 0 such that 

/ u P 2 (-,t) < C(fj,,x,€,P2, IIK, s 0 )||rv 2 > p 2+ 1 (n)) 


for all t G (so,T max ). Thus (|3.3I) is obtained. 


□ 
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An immediate consequence of Lemma 13. II is that Vu is bounded with respect to the norm in L°°(Q). 
Lemma 3.4. Under the assumptions of Lemma \3.‘A for some so G (0,T max ), there is C > 0 such that 


||Vu(-,t)||z,» (n) < C for all t G (s 0 5 ^max ). (3.20) 

Proof. We note that the standard estimate for Neumann semigroup provides ci > 0 such that 

l|Vw(-,t)||x,oo (n) < ||V e t ( A - 1 )u(- jS 0 )IU=o (n) + f* \\Ve^ A -Vu(.,s)\\ L ~ {n) 

J SO 

< e _i || Vu(-, s 0 )||loo(q) + f Ci(1 + {t- s)" = ^^)e _(t_s) ||u(-,s)|| L p 2 ( n)ds 

J S 0 

for all t G (0,T max ). In light of Lemma HTTl we see the conclusion is established. □ 

Now we are in a position to prove boundednes of u. Since the estimate of Vro is nonlocal in time, the 
Moser iteration procedure is still necessary. Upon which, we can finally show that u is bounded. 

Lemma 3.5. Under the assumptions of Lemma \S.2\ for some so G (0,T max ), there is C > 0 such that 

< G f° T all t G (s 0 ,T max ). (3.21) 

Proof. We first see Lemma 13.31 and Lemma 13.41 imply the existence of ci > 0 such that 

lk(-) i )lli=°(n) + l|Vi>(-,f)|| L oo( 0 ) < ci (3.22) 

for all t G (so, Umax)- Testing the first equation in (11.41) with u p ~ l (p > 1), using (12.61) . (13.201) . (13.181) and 
Young’s inequality, we can find constants C 2 ,C 3 > 0 such that 

-4 [ u p + (p- l) f w p-2 |Vu| 2 
P at J n J n 

= (p— l)x [ u p ~ 1 S7u ■ Vu + {p — 1)£ f • X7w + p f u p — p f u p+1 — p f u p w 

< (p — l)x f « p_ 1 Vu-Vr + c 2 ^ [ u p + c f u p v + C 2 f f u p ~ 1 \S7u\+ p f u p 

< P 1 [ u p ~ 2 \Vu\ 2 + (p - l)y 2 f u p \S7v\ 2 + [ u p + c 1 c 2 f [ u p 

4 J n Jn Jn Jn 

+ P [ u p ~ 2 \Xlu\ 2 + c\f 2 p [ u p + p [ u p 

4 Jn Jn Jn 

[ u p ~ 2 \Vu\ 2 + c 3P [ u p 
2 Jn Jn 

for all t G (sojTmax). An obvious rearrangement implies 


d 

dt 



IVlt2 | 2 < C5P 2 



(3.23) 


for all t G (so,T max ) where c 4 ,C 5 > 0 are independent of p. Next, we use (13.231) to perform the classical 
Moser iteration procedure [T] to obtain the boundedness of u. 

Let pk = 2 fe , k G N and M k := sup J n u Pk (•, t). Since Pk > 1, it is easy to find c 6 > 0 such that 

£E(S0 4max) 


d_ 

dt 


,,Pk 



C 4 


\Vu^\ 2 <c 5Pk 2 



u Pk < c e p k 2 / U‘ 


Pk 


(3.24) 
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for all t £ (so,T max ). By means of the Gargliardo-Nirenberg inequality, we get that 


f u Pk — \\u 2 ||| 2 (Q) < C 7 ||Vm 2 ||| 2 ( n )|| , u 2 IlLqfi)) + c t\\ u 2 lllqn) 

J Q 


with a = £ (0,1) and some C 7 > 0 independent of k. Young’s inequality and the definition of p k 

1 i 2 

ensure that there are c§ > 0 and b > 0 satisfying 


c 6 pl / u Pk <c 4 / || 2 + C 8 (pi ) 1 “ 
J O J Q 

<c 4 [ | Vu^] 2 +b k M%_ 1 . 

J n 


C 6 C 7 pl 


ils-')' 


(3.25) 


Combining 


we find that 


4- [ u Pk + [ u Pk < b k Mi 1 . 

dt J a Jn - 

for all t G (sojllnax). The comparison theorem for the above ODE yields 

M k < max{6 fc M|_ 1 , f u Pk (-,s 0 )}. 

Jn 

If b k Ml_ l < J n u Pk (-,So) is valid for infinite k , (13.51) is already derived. Otherwise, a direct induction 
entails 


M k < b k Ml_ x < l2J{k - j) M$ k . 

Taking 2 fe -th root on both sides leads to the assertion. □ 

Now we are ready to prove Theorem [T) 

Proof of Theorem |7J First we see that the boundedness of u and v follow from Lemma 1531 Lemma 13.31 
and (12.31) . Thereupon the assertion of Theorem Q] is immediately obtained from Lemma 12.11 □ 


4 Acknowledgement 

I would like to thank two anonymous referees for carefully reading and giving many suggestions to improve 
the manuscript. 


References 

[1] N. D. Alikakos, L p bounds for solutions of reaction-diffusion equations, Comm. Partial Differential Equations 
4, 827-868 (1979). 

[2] N. Bellomo, A. Bellouquid, Y. Tao, M. Winkler, Toward a Mathematical Theory of Keller-Segel Models of 
Pattern Formation in Biological Tissues, Math. Models Methods Appl. Sci. 25, (2015). 

[3] M. Chaplain, G. Lolas, Mathematical modelling of cancer invasion of tissue: dynamic hetergeneity. Net. 
Hetero. Med 1, 399-439 (2006). 

[4] M. Hieber, J. Priiss, Heat kernels and maximal l p — l q estimate for parabolic evolution equations, Comm. 
Partial Differential Equations 22, 1647-1669 (1997). 

[5] T. Hillen, K.J. Painter, A user’s guide to PDE models in a chemotaxis, J. Math. Biology 58, 183-217 (2009). 


11 




[6] D. Horstmann, From 1970 until present: The Keller-Segel model in chemotaxis and its consequences I, Jber. 
DMV 105 (3), 103-165 (2003). 

[7] D. Horstmann, M. Winkler, Boundedness vs. blow-up in a chemotaxis system. J. Differential Equations 215, 
52-107 (2005). 

[8] K. Osaki, T. Tsujikawa, A. Yagi, M. Mimura, Exponential attractor for a chemotaxis-growth system of 
equations. Nonlinear Anal. 51, 119-144 (2002). 

[9] E.F. Keller and L.A. Segel, Initiation of slime mold aggregation viewed as an instability. J. Theor. Biol., 26 
399-415 (1970). 

[10] T. Nagai, T. Senba, K. Yoshida, Application of the Trudinger-Moser inequality to a parabolic system of 
chemotaxis, Funkcial. Ekvac. 3, 411-433, (1997). 

[11] C. Stinner, C. Surulescu, M. Winkler, Global weak solutions in a PDE-ODE system modeling multiscale 
cancer cell invasion, SIAM J. Math. Anal. 46, 1969-2007, (2014). 

[12] Y. Tao, Global existence for a haptotaxis model of cancer invasion with tissue remodeling. Nonlinear Anal. 
12(1), 418-435 (2011). 

[13] Y. Tao, Boundedness in a two-dimensional chemotaxis-haptotaxis system. arXiv:1407.7382 (2014). 

[14] Y. Tao, M, Wang, Global solution for a chemotactic-haptotactic model of cancer invasion. Nonlinearity 21, 
2221-2238 (2008). 

[15] Y. Tao, M. Winkler, Dominance of chemotaxis in a chemotaxis-haptotaxis model. Nonlinearity 27, 1225-1239 
(2014). 

[16] Y. Tao, M. Winkler, Energy-type estimates and global solvability in a two-dimensional chemotaxis-haptotaxis 
model with remodeling of non-diffusible attractant, J. Differential Equations 257, 784-815 (2014). 

[17] Y. Tao, M. Winkler, Boundedness and stabilization in a multi-dimensional chemotaxis-haptotaxis model, Pro. 
Roy. Soc. Edingburgh, Sect. A 144, 1067-1087 (2014). 

[18] J.I. Tello, M. Winkler, A chemotaxis system with logistic source. Commun. Partial Differential Equations 32, 
849-877 (2007). 

[19] M. Winkler, Aggregation vs. global diffusive behavior in the higher-dimensional Keller-Segel model. .J. Differ¬ 
ential Equations, 248(12):2889-2905, (2010). 

[20] M. Winkler, Boundedness in the higher-dimensional parabolic-parabolic chemotaxis system with logistic 
source. Comm. Partial Differential Equations 35, 1516-1537 (2010). 

[21] M. Winkler, Finite-time blow-up in the higher-dimensional parabolic-parabolic Keller-Segel system, J. Math. 
Pures Appl. 100, 748-767 (2013). 


12 


